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ABSTRACT 


A series system is simulated to obtain lower confidence 
limits on system reliability using Bayesian techniques. A 
comparison between classical and Bayesian methods is made. 
Random beta variate generators are developed and used in the 
simulation. The results of the simulation are tabulated for 
easy comparison of the Bayesian and classical methods. "The 
values of lower confidence limits that are realized using 
the Bayesian method decrease as the number of components 
increase. In most cases, as the number of components 
increase, the Bayesian method appears to yield lower values 


Of lower confidence limits than the classical method. 
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I. INTRODUCTION 


Inas ries system the reliability of the system, RS, 
depends upon the reliability of each of the k components or 
subsystems of which it is composed. This can be modeled as 


follows: 


where Ri 1 με τοι που Of the ith component. Using 
success/fail test data classical or Bayesian methods can be 
ο πο compute Lower confidence limits on RS, In a simple 
system with few components the Bayesian method offers many 
E 5 over the more conservative classical method. 
However, as the system becomes more complex with a higher 
number of components, problems start to arise with the use 
of standard Bayesian techniques. 

Mne purpose of this paper is to compare classical and 
Bayesian methods of computing lower confidence limits.  Com- 
puter simulations were used to determine lower confidence 
limits on system reliability using Bayesian methods. A 
Poisson approximation was used to obtain lower confidence 
Jimits for the classical approach. Several factors were 
porred but the main interest of the investigation was the 
effect of the number of components of a system on the values 


Orelowen econtidence limits on system reliability. 





TTT METHODS 


In order to compute lower confidence limits on system 
liability Of a series system two methods can be used. 
The first method to be discussed is the classical approach 
which is based solely on test results. The second method 
is the Bayesian approach which utilizes assumptions based 


on prior knowledge of similar systems. 


CLASSICAL METHOD 

The classical approach to system reliability uses only 
the mission test data on each component. These results are 
then combined in order to obtain an estimate of the system 
u liability. In the most simple case where k components are 
Exch tested n times and there are no failures in any of the 
Gomponents it is assumed that this is equivalent to testing 
me entire system n times with no failures. In this case 
the lower 100(1-a)% confidence limit can be found in the 
same manner as with just one item. This is done by solving 
lup n the equation: 


Rap = a 


IMS 


MES 


where s = number of successes. 
The solution is a lower 100(1-a)$ confidence limit on RS. 


N 


ma he case of zero failures R5; (a) - ya . When only one 


failure occurs the same procedure can be used. However, 





then more than one failure occurs the procedure becomes 
more complicated. 

In the case where each component is tested n times and 
there are few component failures, an approximation to the 


classical value of RST (o) can be found as follows. 


Let Q,=1-R, , 


~_ 


then RS 


II 


Σ (1-0.). 
DEM "à 


IT the Qi's are small then 


¡St - = the number of failures on the ith component. 
Since fi ic a sum of m Bernoulli tridls 11 is binomial 
(nQ; If Q; is small then each f; is approximately 


Poisson (πο). 


Let F = 


k 
then F is approximately Poisson (n X 0: 
1-1 


mace the upper 100(l-0)% confidence limit [Ref. 1] for 
2 
k η 


ΣΩ. 15 


, where ue 2(F+L) is the 100(1-a)th 
1-1 


2n 


percentile of the x“ distribution with 2(F+1) degrees of 


freedom. 





X2. 2(F*l) 
iaa roximation to the classical value for a lower 


MOOT ge) confidence limit on RS. 


B. BAYESIAN METHOD 
The Bayesian approach to system reliability is to treat 


c i l component reliability, R as a random variable. One 


i? 
particular method is to assume that each component relia- 
bility has a peta distribution. The beta density function 


can be defined as follows: 


F(r:a,b) = ETET Yarleı-ny Pl for 0 <r <1, a >0, b > 0, 


where B(a,b) = ο, 


prior beta density B(r;a,b) is assumed for E component 
before testing begins. The test data is then used to obtain 
ED osterror beta density B(r;a*s,b*f) where s is the number 
Of successful tests and f is the number of failures for the 
component. 

By generating a random beta variate for each component 
based on its posterior density a value of RS is calculated 


using the model: 


A distribution of RS values is simulated by repeating this 


procedure. The lower confidence limits are thus estimated 





by the appropriate percentile points of the simulated 
dutributions of RS values. 

The main advantage of the Bayesian method is that past 
experience with technically similar hardware can be used to 
determine appropriate priors for each component. This allows 
Meher reliability goals to be met with less testing. How- 
Per, Some precautions must be observed in selecting the 
beta priors. Choosing a beta prior of B(r;a,b), where a and 
b are integers, is equivalent in the classical sense to 
assuming that (a+b-1) tests have been performed and that a 
of these were successes, It is equivalent because the ath 
percentile point of the B(r;a,b) distribution is the lower 
100(1-a)$ confidence limit for p when there are (a*b-1) 
Bernoulli trials, each with probability p of success, and a 
of these trials are successful. This is readily apparent 
from an identity that relates the beta distribution to the 
ermomial distribution [Ref. 2]. 

τ ΠΤ; ΠΕ Dela parametels (atb) in the prior 
density increases relative to the number of tests to be 
observed the resulting lower confidence limits are more of 
Æ unction of the assumed prior density than of the test 
mestilts, One procedure is to limit the sum of the two 
Parameters to a function of the number of tests. If the b 
Em meter is less than one a different problem exists even 
though the sum of the parameters is small. This problem 


can be shown by the following discussion. 





ἵ η ο ὃ Ρος peta distribüutions exist such 
ο πο ο crcentalesposnbort each 1s the same. If the 
b parameter is set equal to 1, then for any B(r;a,b) 
ο του an a* can be found such that the ath percen- 
NES cu OL B(r;a,b) and B(r;a*,1) are the same. The 
B(r;a*,l) ο ο bution Is convenient since it can be related 
the classical case of having already tested a^ items with 
pemrsrlures. The problem when the b parameter is less than 
one can be shown by converting the initial B(r;a,b) prior to 
T31 1l) prior, In this case if there are no failures 
mie posterior distribution will be B(r;a*+s,1) where s is 
|| number of successful tests. If the initial beta prior 
15 not converted the posterior distribution will be 
B(r;a+s,b), Koach can benconverted to B(r; (ats) ο). In 
all cases that were investigated the (a*s)* was larger than 
dus when the b parameter was less than one, Thus, in an 
divalent classical sense, it appears that each successful 
mre will be counted as more than one success. This 
amplification effect seems to increase rapidly as the b 
parameter approaches zero. Also if the beta prior (B(r;a,b), 
where b is less than one, is converted to a beta B(r;a*,1) 
momcmisure that eachysuccess counts only as ome success, the 
fesulting a* may be too large. The details of the conversion 
are given in the next section and a table of the resulting 


us for various values of b can be found in Appendix A. 





III. SIMULATION 


The purpose of this investigation was to compare lower 
confidence limits derived by using Bayesian methods with 
those obtained by using the classical method. In order to 
simplify the simulation and to be able to make comparisons 
it was assumed that each of the components had the same 
beta prior density and that each component was mission 
tested the same number of times. Using the same beta prior 
for each component is not a necessary condition for the 


simulation to work, 


A. SIMULATION PROCEDURE 

A series system of k components was simulated in the 
following manner. Let Bl.(r;a,,b.) be the initial beta 
ior for R,. For this investigation the following three 
mara priors were used:  B(r;5.0,0.05), B(r;5.0,0.05) and 
DES, 5.0,0.10). For each Bl, (r3a;,6;) a new beta prior, 
pr ,a,1), was Computed Such that the l0th percentile 
Paint of ΗΓ was the same as the 10th percentile 
meant of Bli(r;ai,bi). The value for a: is easily computed 
mece if the b parameter is 1 then the percentile point (P) 
HFA B(x:a, l) distribution is given by the equation P = xo. 
Thus a; ANA where x 1s theiloth percentile 


ponunt of Bli(r;a;,b:). 


tl 


LEE n.: chesmumbersotstrials for cach component; 


1 


and Si the number of successes for each component. 


10 





Tne B2. 1{1; a^ Ec x πει adjusted in the following 


manner: 
ba 1{1; 81 1) where QE min(0.75n; ,a2) 
B22,(r;a;,1) where acis min(1.0 || να.) 
B23,(r;a;,1) where a, 7 mias lys) 
A beta posterior B5; (r;ai,bi) for each R; was computed in 


Tour cases as follows: 


Case 1 B3. )نت‎ a; De 1) 


BMT aT botes) 
1 ΠΠ =] 


Case 2 B3. j (r; a; lo): i) 


B21. ن‎ ats. ,l*n:i-S; ) 


Case 3 B3. ία; a: PU i) 


BI2. ; C a ,l*ni-S; ) 


Case 4 B3. ; (7; a; "Da i? 


Do Ure cM sS 
i JS J 1 


By generating B3. quee a; abe i) random variates, which represent 
the posterior distributions of each of the components, a 
pue for RS was obtained by taking their product. This 
pusccdure was repeated 500 times and realized values for RS 
were then ordered. The lower 100(1-a)% confidence limits 
.', Den determined by selecting the appropriate percentile 


points of the simulated distribution on RS. 


B. GENERATION OF RANDOM BETA VARIATES 

Two means of generating random beta variates were used 
in the simulation. In the case where the b parameter of the 
B(r;a,b) distribution was an integer an exponential genera- 


OT vas Used tO realize random beta variates. In the case 


11 





of noninteger parameters, Monte Carlo rejection techniques 
were used. 

In order to generate random B(r;a,b) variates where b 
is an integer the following logic was used [Ref. 3]. Assume 
ES B(v;a,b). 
Define U - -1nY. 


Then the moment generating function 


M, t= Ele’ ] = Ele” ]= EN η 
Thus 
1 
ران‎ 1 a-t-1 0-1 
POENI E B(a,5) y » AO y 


Blas tb) τα afb) l (a=) 
p a, . F(3)-T(5) as host 





tl 

3 

7 

t| — 
wirt 

os 


Since c=) 
D: 


is the moment generating function for an exponential random 
variable with parameter A, then U is the sum of b indepen- 
dent exponential random variables with failure rates a, 


ERU وهن‎ Ab, 


MC. SE Us where U: 1s exp (a+j). 
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πες Y = E random beta variates can be generated with 


the use of existing exponential generators. For the simu- 
lation, where the posterior distributions on R; were 
B3.(r;a;,b.) απο ۴ E were integers, values for RS 
were obtained by using the following formula: 
k b.-1 
-]nRS = 2 ; دول‎ : 
i=l j=0 J 


where 


E is exponential (a;+j). 


In the case where the b parameter of the beta density 
was noninteger, random beta variates were generated by 
mame Monte Carlo rejection techniques [Ref. 4]. Since 
this investigation was concerned with small ee ος thE SD 
Parameter the Monte Carlo technique had to be modified 


^ 


Slightly. The procedure that was used is outlined below. 


1. Compute P | ١ 11601 < OO) 01676 το δα ο. 
2. Set H H=Max value of fx (x) Lone s acme ο 
3. Generate R R=Random U(0,1) number. 

Ií R > P Then X.=1.0, Go to 3. 

s Generate RI | RI-Random U(0,1) number. 

ο. If R1 > 0.999 μου,  οο τσ 

fee Generate R2 R2=Random U(0,1) number. 


ο Li (RZ < fy (R1) Then X,=Rl, CO UO 


9. Go to 5. 


This procedure was repeated until 500 random beta variates 


were realized for each posterior density. 
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C. MODIFICATION OF SUBROUTINE BDTR 

Subroutine BDTR [Ref. 5] was used in the simulation to 
convert beta priors and to compute beta densities in the 
generation of random beta variates using Monte Carlo tech- 
Dates. The subroutine BDTR computes the probability that 
the random variable U is less than or equal to x, where U 
is distributed according to the beta (u;a,b) distribution 
ο «χε |. The value of the density at x is also com- 
puted. In order for the computations to be valid the sum 
of the parameters must be greater than or equal to one. 
Mm the subroutine BDIR this condition is met by restricting 
botn parameters to values greater than or equal to 0.5, 
The subroutine was modified so that smaller values of the 
b parameter could be used since the a parameter was always 
@meater than or equal to one. Two additional variables 
were added to the parameter list of the subroutine which 
miemeased the efficiency of the simulation by acting as 
flags and saving values when the subroutine was being 
called during the generation of random beta variates. All 
modifications to the subroutine are shown in the listing 


of the computer program. 
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IV. RESULTS 


The results of the simulation are given in tables for 
ο ο ος the three beta priors that were used. The number 
of failures was the total number of failures in all of the 
components. For n=50, and in a similar manner where n=75 


and n=100, the failures were assigned as follows: 


Failures Successes 

0 5, 3 SO اوو مول ول هم‎ ٢ 
1 S; = 49, i=l 

SO MESE UC NA 
2 SENS ar د‎ 

= 50, 155,4, SK 
3 ج‎ 48, i=l 

= 49, i=2 

πο as K‏ اه 
S 179‏ 4 

= 50, 1=3,4, ,K 
5 In A س‎ 

= 49, 153 

2550573 sd 5, v vd 
6 ΠΕ 45 2 

= 49, 155,4 

xU Eq OS ue 


The beta priors for the Bayesian cases were as follovs: 


Case 1 Bl(r;a,b) 


Case 2 Brisa.) a πιπ 7 sna) -n=50 775.100 


Case د‎ ۹22) 3,1 3 19۱11) .0 H.A) N 50,15,100 


Case 4 B25 (rsa bea 


πη 0 
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Me appropriate values for Bl(r;a,b) and BZ(r;a*,1) are 


given at the top of each table. 
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90% 


Bl(r;5.0,0.03) 


Lower Confidence Limits on RS 


converts to 


Fifty tests for each component 


k Failures 


40 
40 
40 
40 
40 
40 
40 
30 
30 
30 
30 
30 
30 
30 
20 
20 
20 
20 
20 
20 


20 


0 


Classical 


Case 

.9540 
29222 
Bo 
. 8664 
. 8401 
.8145 
.7894 
. 9540 
.9222 
20095 
.8664 
.8401 
‚8145 
.7894 
.9540 
79222 
.8936 
. 8664 
. 8401 
.8145 


22595 


Case 1 
0.9532 
0 ٢ τ2 
0.9038 
0.8827 
0.8647 
0.8410 
0.8247 
0.9616 
0.9345 
0.9119 
0.8897 
0.8696 
0.8491 
0 د‎ 7 
009751 
0.9398 
0.9159 
0.8950 
0.8754 
0.8549 


0.8548 
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B2(r;596.8,1) 


Bayesian 
Case 2 Case 5 
EI 016162 
5759 0 6151 
eo) ΠΟ; 
‚5541 0.5966 
.5493 0.5920 
. 5429 0.5861 
. 55358 0. 5775 
. 6547 ο ο ον 
.6467 060829 
.6384 0.0752 
02968 0.6636 
6247 020627 
. 6179 03 
. 5 0.6467 
HRS 2 0" 
2/3/00 01035 
727 ο τοι 
1022 0.7495 
.7074 0.7388 
. 7012 ) د‎ ٢ 
. 6924 0 7251 


Case 4 
0.6789 
0.6779 
0.6663 
0.6116 
U 06575 
ο ο. 
0.6446 
0.7434 
057571 
0.7305 
0.7205 
0.71٨٢ 
0.7140 
0 7057 
0.8139 
0.8074 
0.8006 
0.7940 
0. 7850 
077 


ο τος 





902 Lower Confidence Limits on RS 


ET S0005) converts to —B2(r;596.8,1) 


Fifty tests for each component 


Classical Bayesian 
k Failures Case Case 1 Case 2 Case 3 Case 4 
10 0 0.9540 . 9840 . 8478 0.8654 .8908 
10 1 050222 . 9474 ‚8413 0.8597 .8861 
10 2 0.8936 qo . 63506 0.8501 πολ 
10 5 0.8664 .9023 .51506 0.8368 2510 725 
10 4 0.8401 6617 .8069 0.8289 . 8607 
10 S 0.8145 ο. © .7949 0n SIS . 9538 
10 6 0.7894 . 8399 207 7 0.8067 .8423 
5 0 0,1 “9924 ο ο 0.9264 . 9407 
5 t ο ο” nop 29075 0.9140 29506 
5 2 DI . 90259 . 8850 0.8986 ο. 
5 3 0.8664 . 9062 8105 0.8910 29d 19 
5 4 0.8401 .9833 . 00657 08738 . 9027 
5 5 0.8145 .8616 . 8466 0.8645 ου; 
5 6 0.7894 .8435 .8400 0.8587 . 8856 


 w——— rs am rr Fw — s T — < -— —— — E — — = — —— — m 
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902 Lower Confidence Limits on RS 


νυ η) converts to  B2(r;596.8,1) 


Seventy-five tests for each component 


Classical Bayesian 

k Failures Case Case 1 Cas 2 Case š Case 4 
40 0 N ο. ται 0.6915 . 7242 ΙΤ 
40 1 0.9481 .9496 0.6906 222553 0. 7717 
40 2 υπ" 2956505 0.6794 7/150 1029 
40 038209 nons 0.6747 . 7088 7593 
40 4 0.8934 0052 0.6708 “2051 25551 
40 5 013 18 265 0.6658 .7006 20573 
40 6 0.8596 5721 0"6552 ο ο . 7463 
50 0 ο το” 207756 0.7540 ne .8206 
50 1 0.9481 299 0.7478 2 5 . 8160 
30 2 0.9290 . 2 0.7415 HO -STI 
30 5 0 0105 .9186 0. 7315 . 7610 . 8038 
30 4 0.8934 2900 2 00 2200.2 . 6031 
30 5 0.8763 .8920 077256 -7553 τος 
50 6 085906 .8780 TS .7480 240287 
20 0 021098 .9803 Ον ΘΠ .8425 ο le 
20 1 0.9481 . 9591 0.8157 -8367 .8671 
20 2 0.9290 0455 0 5091 . 8309 8022 
20 3 0" 105 0235 0.8028 2025] σας 
20 4 0.8934 .9093 0.7941 .8174 .8510 
20 5 0.8763 03958 0.7895 noy .8476 
20 6 νυ ΝΠ .8804 0.7829 5 ποίη 
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90% Lower Confidence Limits on RS 


ιν. 5 0 0 05) converts to B2(r;596.8,1) 


Seventy-five tests for each component 


Classical Bayesian 
k Failures Case Case 1 Case 2 Case 3 Case 4 
10 0 0.9693 0.9891 0.8957 0.9082 0.9258 
10 1 .9481 .9647 0.6912 . 9041 ΟΙ 
10 2 292 90 .9474 0.8836 . 8974 οι τι 
10 3 SONS 5 0 67352 .8882 2095 
10 4 .8934 .9138 0.8669 0925 .9049 
10 5 .8763 9001 0.8584 2507250 8987 
10 6 ‚8596 .8836 0.8492 .8668 . 8920 
5 0 . 9693 .9940 0.9434 . 9503 . 9600 
5 1 . 9481 OI 0.9338 .9418 .9532 
5 2 . 9290 . 9504 0.9218 . 9515 20m 07 
5 S .9109 :9:3:04 0 91060 . 9261 .9405 
5 4 oe oe 9 0.9072 ‚9183 . 9342 
5 5 ο ος 29002 0.8951 . 9077 3255 
5 6 ‚8596 .8850 0.8908 .9038 .9224 
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90% Lower Confidence Limits on RS 


1 s 050,035) converts to B2(r;596.8,1) 


One hundred tests for each component 


Classical Bayesian 

k Failures Case Case 1 Case 2 Case 3 Case 4 
40 0 0.9770 MO 27:9 9:5 0270850 3299 
40 1 0.9611 .9586 5 ο ος 507 51 
40 2 0.9458 . 9481 .7483 0.7759 81.65 
40 3 0 5552 9560 . 7445 0 7725 154 
40 4 0.9200 .9266 ΠΠ 0.7695 .8109 
40 5 0.9073 ‚9106 20772 0 5 . 8079 
40 6 0.8947 ΠΡ 15083 0.7601 .8030 
30 0 09770 57 .8091 028508 28022 
30 1 ΤΙ EOD .8042 0.8264 05095 
30 2 0.9458 ¿9514 BE 0.8218 Sg 
30 5 0.909552 .9417 .7913 0.8149 .8490 
30 4 0.9200 29270 . 7906 ο ο .8484 
30 5 0.9075 4 تنا‎ .7862 0751 7 ΙΙ 
50 6 0" 1; -9070 ΠῚ 0.8043 EU 
20 0 ο 770 οσο . 6052 0.8793 ο 
20 1 0.9611 .9685 55295 0 5749 . 8986 
21 2 0.9458 οτι . 9552 0.8703 .8948 
20 3 ο. ο .9446 .8481 0.9658 . 8911 
20 4 0.9200 η nodus 0.8596 .8861 
20 5 0:9 0/755 9165 89557 0.8564 . 4 
20 6 0.8947 . 9090 8324 0.8518 257096 


41 





909 Lower Confidence Limits on RS 


INE 5 020.05) converts to  B2(r;596.8,1) 


One hundred tests for each component 


Classical Bayesian 
k Baa lures Case Case 1 Case 2 Case 3 Case 4 
10 0 0.9770 ‚9918 ο ο 0.9303 . 9438 
10 1 0.9611 29027 0.9772 0 9272 .9415 
10 2 0.9458 . 9595 0 1 0.9220 SAL 
10 5 0.9332 .9474 0023054 0.9150 49:514 
10 4 υ ου 29550 0.8984 0.9106 92718 
10 5 0.9073 .9210 0.8919 0.9047 2:90:27 
10 6 0.8947 AOT 0.8847 0.8984 75 
5 0 0:9770 -9956 Omen 17 029025 „9699 
5 1 υπ οι -9753 0.9499 0.9560 . 9647 
5 2 0,9458 -9009 0.9408 0.9480 . 9582 
5 3 ο ο” .9504 0.9364 0.9441 .9551 
5 4 079200 . 9344 0.9296 0.9582 2:9 505 
5 5 0 5075 92702 0.9204 0.9300 .9436 
5 6 0.8947 . 9121 0:9171 0 0271 .9413 


τε 





90$ Lower Confidence Limits on RS 


ET 5.0005) converts to .B2(r;136,8,1) 


Prey tests for each component 


Classical Bayesian 

k Failures Case Case 1 Case 2 Case 3 Case 4 
40 0 0.9540 9344 SOL 0.6162 .6789 
40 1 0 0222 ΙΙ» 25039 0.6٢ Soa 
40 2 USO .9844 9399 0.6020 .6663 
40 S 0.8664 η ‚5541 0.5966 .6116 
40 4 0.8401 . 8459 .5493 05320 .6575 
40 5 0.8145 .8265 . 5429 0.5861 0525 
40 6 0.7894 .8091 5538 SS .6446 
30 0 0.9540 .9441 6547 0.6903 ~ 7434 
30 1 DET . 9203 . 6467 0.6329 SS 
30 2 0.8936 ‚8918 .6384 0. 6752 . 7305 
30 3 0.8664 . 9758 .6258 0.6636 ο» 
30 4 0.8401 . 8542 ‚6247 0.6627 . 7196 
30 5 0.8145 SONS . 0179 0.53 .7140 
30 6 0.7894 0155 .6075 0.6467 .7057 
20 0 0.9540 . 9608 .7452 0. 7731 61595 
20 1 079222 2η . 7366 ο του. .8074 
20 2 01.77۴ .9072 MULTI ο 5.7 . 8006 
20 3 0.8664 .8874 . 7192 0.7495 . 7940 
20 4 0.8401 . 8607 .7074 0.7388 .7850 
20 5 0.8145 .8385 SU EZ ο ο. .7802 
20 6 0.7894 26213 .6924 097251 .7734 
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SU tower conridence Limits on RS 


Hu 00205) “converts to .B2(r;136.8,1) 


Fifty tests for each component 


Classical Bayesian 
k Failures Case Case 1 Case 2 Case 3 Case 4 
10 0 0.9540 0 47 0.8478 0.8654 0.8908 
10 1 9222 0.9452 0.8413 0.8597 . 8861 
10 7 8350 I? 0.8306 0.8501 . 8782 
10 3 . 8664 .8934 0.8156 0.8368 .8673 
10 4 . 8401 802 0.8069 0.8289 .8607 
10 9 .8145 .8469 0.7949 0.8182 2555/6 
10 6 .7894 -8313 077872 0.8067 ‚8423 
5 0 . 9540 . 9866 0.9163 0.9264 ΟΠ 
5 1 02 non 0 ۹243 0.9140 οσο 
5 2 .8936 29254 0.8850 0.8986 . 9180 
5 3 . 8664 ΠῚ 028703 0.8910 . 9119 
5 4 .8401 28125 057 0 37 290 7 
5 6 694 . 84 0.8400 0.8587 . 8856 
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902 Lower Confidence Limits on RS 


NT 5.0,0.05) converts to B2(r;136,8,1) 


Seventy-five tests for each component 


Classical Bayesian 

k Failures Case Case 1 Case 2 Case 3 Case 4 
40 0 09695 553 0 35 υ د‎ . 1 
40 1 . 9481 388 0.6906 0 7255 a7 
40 2 .9290 50209 0.6794 057150 7023 
40 3 . 9109 .9046 0.6747 0.7088 159395 
40 4 .8934 .8943 0.6708 0.7051 ο 
40 5 ο .8769 0.6658 0.7006 “1525 
40 6 ‚8596 .8601 0 072 0.6936 .7463 
30 0 .9693 .9615 0.7540 0.7811 . 8206 
30 l .9481 .9449 0.7478 027755 .8160 
30 7 .9290 9310 0. 7415 0.7697 poH TE 
30 3 .9109 9099 0 537 027610 . 8038 
30 4 uo .8988 0 7510 0.7602 8051 
30 5 .8763 . 8808 022256 0. 7555 . 7989 
30 6 . 8596 .8684 0.7175 0.7480 ο. 
20 0 .9693 .9710 028219 0.8423 δα 
20 1 .9481 89526 ΟΡ 065067 8071 
20 2 .9290 .9367 0.8091 0.8309 „8622 
20 3 . 9109 20172 0.8028 525 10.575 
20 4 . 8934 "0057 0.7941 0.8174 δυο 
20 5 SAU .8885 0.7895 0.8132 .8476 
20 6 .8596 25770 0.7829 0.8073 . 8427 
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127 5.0,0.05) 


90% 


016۳٣٤٣ LO 


B2(r;136.8,1) 


Seventy-five tests for each component 


k 
10 
10 
1 
10 
10 
10 
10 


Failures 
0 
1 


Classical 


Case 
0.9693 
0.9481 
0.9290 
0.9109 
0.8934 
0.8763 
0.8596 
0. 9693 
0.9481 
0 0 050 
.9109 
. 8934 
0.8763 


02330 


Case 1 
0.9826 
0 5597 
0.9421 
0.9245 
0.9101 
0,8959 
0.8819 
0.9892 
0.9657 
029476 
0.9280 
0 31 
0.8972 


0.8842 
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Bayesian 


Case 2 
07 
0 5912 
0.8836 
O. 6752 
0.8669 
0.8584 
0.8492 
0.9431 
0.9338 
0.9218 
0.9160 
050072 
0.569551 


0.8908 


Lower Confidence Limits on RS 


Case 3 
0.9082 
0.9041 
0.8974 
0.8882 
ο ο 
0.8750 
0.8668 
0.9503 
0.9418 
O. 
015261 
0.9 3 
023077 


020055 


Case 4 
0.9258 
0 9225 
0.9171 
0.9095 
0.049 
0 377 
0.8920 
0.9600 
Ὁ 
0.9447 
0.9405 
0.9542 
0 5 


0.9224 





902 Lower Confidence Limits on RS 


ιν, 50.0.05) converts to B2(r;136.8,1) 


One hundred tests for each component 


Classical Bayesian 
k Failures Case Case 1 Case 2 Case 3 Case 4 
40 0 ου 0.9649 0.3 0.7850 0.8151 
40 1 . 9611 0 39507 υ ο. 0.7843 Π.Ι; 
40 Z . 9458 093573 0.7483 0.7759 0.8071 
40 Š 29597 00277 0.7445 ο 555 0.8015 
40 4 SOM) 0.9148 0.7412 μυ. ο 0.8015 
40 5 ο. 0.9086 07272 0272058 0.7985 
40 6 .8947 0.8945 0.7308 0.1 O70 
30 0 . 9770 029721 0.8091 0.8591 0.8308 
30 1 .9611 0.9579 0.8042 0.8264 0785135 
30 2 .9458 0.9444 O 7991 0.8218 0.8473 
30 B 29557 0.9336 0.7913 0.8149 0.8413 
30 4 . 9200 0 307 0.7906 0.8142 0.8400 
30 5 .9073 0:01 0 0.7862 0.8102 0.8372 
30 6 .8947 0.9011 T 0.8043 0.383515 
20 0 τι IEA 0.8632 078793 0.8970 
20 1 .9611 0.9641 ο ο ὃν 0.8749 0769352 
20 Z . 9458 0194906 0 2 ο 505 0.8895 
20 5 29052 0.909595 0.8481 0.8658 0.8854 
20 4 . 9200 0.9270 0.8413 0.8596 0.8801 
20 5 290755 0290172 05 377 0.8564 0 4 
20 6 . 8947 0.9063 0 4 0.8518 0598738 
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902 Lower Confidence Limits on RS 


SO 005) converts to  B2(r;136.8,1) 


One hundred tests for each component 


Classical Bayesian 
k Failures Case Case 1 Case 2 Gase 3 Case 4 
10 0 09770 -9871 ο οἱ . 9303 . 9408 
10 1 0.9611 . 9690 ο د‎ 2272 .9382 
10 2 0.9458 .9566 0.9114 3220 293338 
10 3 09352 . 9457 0.9034 . 9150 9277 
10 4 0.9200 9527 0.8984 ΠΕ . 9240 
10 5 059075 . 9219 (S SIE .9047 . 9190 
10 6 0.8947 79091 0.8847 . 8984 9135 
5 0 ΠΠ 2519/95 υ 9025 . 9082 
5 1 0.9611 3759 0.9499 . 9500 9027 
5 2 0.9458 19599 0.9408 .9480 . 9560 
5 3 0.9332 . 9478 0.9364 . 9441 .9526 
5 4 0.9200 29 9 0.9296 .9382 . 9476 
5 5 0.9075 . 9271 0.9204 .9300 . 9406 
5 6 0.8947 314 0 9171 O οσοι 
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902 Lower Confidence Limits on RS 


ο. 5 '0 0110) converts to B2(r:39.5,1) 


Fenty tests for each component 


Classical Bayesian 
k Failures Case Casel Case 2 Case 3 Case 4 
40 0 .9540 0.8825 Forol 0 نن‎ 1 Same 
40 1 50222 0.8639 295459 0.5809 mes 3 
40 2 229530 0.8407 . 5599 0 5071 
40 5 . 8664 0.8230 . 5541 0559155 
40 4 .8401 0.8028 „493 075506 
40 5 .8145 0 7672 . 5429 0.5505 
40 6 0.7894 0.7666 075353 0 5 T 
30 0 0.9540 0.9055 0.6547 0.6608 
30 1 0.9222 0.8842 0.7 0.6529 
30 D .8936 0.8619 σος 0.6447 
30 3 .8664 0.8423 POZOS 00322 
30 4 . 8401 (ST 6247 0.063515 
30 5 .8145 08073 .6179 0.6179 
30 6 . 7894 ο 4955 .6075 0.6142 
20 0 .9540 029515 πας, 0.7500 
20 1 noz 0.9083 .7366 0.7416 
20 2 . 8936 0.8864 ran 0, 73:28 
20 3 . 8664 0.8633 ‚7192 0.7244 
20 4 .8401 0.8414 . 7074 021128 
20 5 .8145 ο οἱ .7012 0. 7067 
20 6 .7894 0.8040 .6924 0.6981 
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90% Lower Confidence Limits on RS 


EE 5 050-10) converts to B2(r;39.5,1) 


Fifty tests for each component 


Classical Bayesian 

k Failures Case Case 1l Case 2 Case 3 Case 4 
10 0 0.9540 0.9600 0.8478 0.8508 Same 
10 1 029222 0.9360 0.8413 0.8445 esc 3 
10 2 79956 .9084 . 8306 0.8340 
10 3 .8664 15592 .8156 0.565155 
10 4 . 8401 . 8645 . 8069 0 , 5 ۰7 
10 5 .8145 .8456 . 73949 0.7990 
10 6 .7894 .8230 NU 0.7865 

5 0 . 9540 SN SOME 0.11 

5 1 ο... . 9446 .9023 0.9044 

5 2 2595/5 9207 ΕΤ 0.8873 

5 3 .8664 SONA τος 0.8741 

5 4 .8401 δε. 35555 / 0.8665 

5 5 .8145 5551 .8466 0.8497 

5 6 .7894 8304 .8400 0.8433 
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90% Lower Confidence Limits on RS 


TUNE. 50.0.10) converts to B2(r;39.5,1) 


Seventy-five tests for each component 


Glassical Bayesian 
k Fallures Case gases Case 2 Case 5 Case 4 
40 0 ου ον 0.9193 026551 Same Same 
as as 

40 1 0.9481 0.9041 0.6541 Case. 2 Case e2 
40 2 . 9290 615 0.6419 
40 3 πο 6755 0.6369 

40 4 . 84 . 8641 0.86 
40 : .9763 29520 0 6 71 
40 6 ΟΞ 28305 0.6190 

30 0 .9695 29577 0.7234 

30 1 . 9481 . 9201 ο. προ 

30 2 . 9290 .9044 0.7096 

30 3 .9101 . 8884 0.6990 

30 4 . 8954 20002 0.6981 

30 5 20165 . 8626 0.6922 

30 6 55 29:509 0.6834 

20 0 S .9541 0.7986 

20 1 . 9481 19:597 0.7916 

20 2 „9290 49221 0.7843 

20 3 ΙΙ SES 057775 

20 4 .8934 .8925 07677 

20 5 "57/05 . 8800 0.7625 

20 6 530 . 8616 0. 7552 


31 





90$ Lower Confidence Limits on RS 


TT ST 00.10) converts to B2(r;39.5,1) 


Seventy-five tests for each component 


Classical Bayesian 
k Failures Case Case 1 Case 2 Case 3 Case 4 
10 0 0.9693 0.9718 0.8814 Same Same 
as as 

10 1 0.9481 0.9520 0.8763 Case 2 Case 2 
10 2 SD 0.9369 0.8677 
10 3 ΠΠ 0.9167 0 5559 
10 4 . 8954 0.9057 0.8488 
10 5 “81703 0 6912 0: 5595 
10 6 . 8596 0.8758 0.8290 

5 0 BS 0 7 0.9354 

5 1 .9481 0.9604 0.9244 

5 2 EIU 0.9446 0.9108 

5 3 τ 079257 0.9042 

5 4 .8934 0.9122 0.8942 

5 5 29503 0.8989 0.8806 

5 6 . 8596 ν᾽ 0.8750 
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90% Lower Confidence Limits on RS 


ιο 5 ϱ 0.10) converts to B2(r;39.5,1) 


One hundred tests for each component 


Classical Bayesian 
k Failures Case Case 1 ο ο” Case 3 Case 4 
40 0 0 0770 0.9411 7067 Same Same 
as as 
40 1 0.9611 0 ۹6 .7058 Case 2 Case 2 
40 2 0.9458 0.9130 .6950 
40 3 0.9332 0.9045 .6905 
40 4 0.9200 0 5525 .6867 
40 5 Ον 0.8855 .6819 
40 6 0.8947 0.8743 non 7 
30 0 0 5770 0.9495 .7666 
30 1 0.9611 (29/377 τους 
50 2 0.9458 09258 22527 
30 3 0,9552 0.9164 es 
30 4 0.9200 0.9034 . 7446 
30 5 0.9073 0.8960 .7394 
30 6 0.8947 0.8838 57 
20 0 0.) 0.9652 46:915 
20 1 0.9611 0.9517 08255 
20 2 0.9458 0.9379 ‚8193 
20 3 079332 0 9276 13597 
20 4 0.9200 0.9175 .8050 
20 5 0.9073 0.9074 .8006 
20 6 0.8947 0.8969 0.7943 
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902 Lower Confidence Limits on RS 


135 -0.0.10) converts to B2(r;39.5,1) 


One hundred tests for each component 


Classical Bayesian 
k Failures Case Case 1 Case 2 Case, Case 4 
10 0 ο ο Ον 0.9016 Sane Same 
as as 

10 1 ΙΙ ΙΙ; 0 35 Case ñ? Case 2 
10 2 (9m 385 1-157 0.8901 

10 3 191597 .9409 0.8802 
10 4 . 9200 29780 0.8742 
10 5 90075 .9148 0.8661 
10 6 .8947 . 9062 0208375 

5 0 ο. τ; 0.9466 

5 1 . 9611 200 0,9570 

5 2 .9458 25550 0.9262 

5 3 29552 . 9450 09207 

5 4 00 acr 0.9124 

5 5 E :9217 0.9010 

5 6 . 8947 ΙΖ 0 59065 
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V. =EONEEUSTONS 


The results of the simulation indicate that the classical 
ο σας) to System reliability usually yields higher values 
of lower confidence limits than the Bayesian approach. When 
there are 10 or more components, only the beta priors with 
b parameters less than one provided higher values of lower 
confidence limits than the classical approach. When using 
meee bavesian approach, the effect of each prior density 
assumption is to obtain a more optimistic posterior density 
ön each component reliability. However, there is some proba- 
bility mass that is still assigned to small values of each 
Ri. Thus, as the number of components increase the probabil- 
"ος small values of RS increases since RS is a product of 


the R;'s. Therefore, as the number of components increase, 


Meee fect Of the prior density assumptions tends to diminish. 
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APPENDIX A 
BETA DENSITY CONVERSIONS 


The following table lists values of a*, where a* is such 
mae the 10th percentile point of a B(x3;a,b) distribution is 
the same as the l0th percentile point of a B(x;a*,1) distri- 
lll ion. 


Table of a*!s 


b a=1 a=5 a=10 

100 00 5.0 10.0 

0.90 دازا‎ 5.506 ος 
0.80 1.10 5 71 JI erm 
0.70 ley 0:597 1 
0.60 1. 4 2209 14.46 
ο 50 DES 83: 10.01 
0.45 1.47 87 142 
0.40 m SO les 
0 5د‎ ToS 10.68 24-7 
0.30 ΠΣ; 17-7 25515 
0:25 10 14.23 23235 
0.20 2,56 77236 50.57 
0.18 202 J 5 40.70 
0.16 3216 ZZ 46.16 
0.15 Soe 25575 49,99 
0.14 Seel 25.509 55.09 
0 17 4.28 30.96 64.78 
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